NMR evidence for a strong modulation of the Bose-Einstein Condensate in BaCuSi20e 



O 

o 

(N 
5h . 

< 



S. Kramer. 1 R. Stern, 2 M. Horvatic, 1 C. Berthier, 1 T. Kimura, 3 and I. R. Fisher 4 

1 Grenoble High Magnetic Field Laboratory (GHMFL) - CNRS, BP 166, 38042 Grenoble Cedex 09, France 
2 National Institute of Chemical Physics and Biophysics, 12618, Tallinn, Estonia 
s Los Alamos National Laboratory, Los Alamos NM 87545, USA 
J 'Geballe Laboratory for Advanced Materials and Department of Applied Physics, Stanford University, Stanford CA 94305, USA 

(Dated: February 5, 2008) 

We present a 63,65 Cu and 29 Si NMR study of the quasi-2D coupled spin 1/2 aimer compound 
BaCuSi2C>6 in the magnetic field range 13-26 T and at temperatures as low as 50 mK. NMR data in 
the gapped phase reveal that below 90 K different intra-dimer exchange couplings and different gaps 
(Ab /Aa = 1.16) exist in every second plane along the c-axis, in addition to a planar incommensurate 
(IC) modulation. 29 Si spectra in the field induced magnetic ordered phase reveal that close to the 
quantum critical point at H c \ = 23.35 T the average boson density n of the Bose-Einstein condensate 
is strongly modulated along the c-axis with a density ratio for every second plane ha/tib — 5. An 
IC modulation of the local density is also present in each plane. This adds new constraints for the 
understanding of the 2D value (f> = 1 of the critical exponent describing the phase boundary. 

PACS numbers: 75.10. Jm,75.40.Cx,75.30.Gw 
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The interest in Bose-Einstein condensation (BEC) has 
been considerably renewed since it was shown to occur 
in cold atomic gases p|. In condensed matter, a formal 
analog of the BEC can also be obtained in antiferromag- 
netic (AF) quantum spin systems 0, d, 0, under an 
applied magnetic field. Many of these systems have a 
collective singlet ground state, separated by an energy 
gap A from a band of triplet excitations. Applying a 
magnetic field (H) lowers the energy of the M z = — 1 
sub-band and leads to a quantum phase transition be- 
tween a gapped iron magnetic phase and a field induced 
magnetic ordered (FIMO) phase at the critical field H c \ 
corresponding to A mnl -g^BH c i = 0, where A m ; n is the 
minimum gap value corresponding to some q vector q m i„ 
@; S EL Q • This phase transition can be described as a 
BEC of hard core bosons for which the held plays the role 
of the chemical potential, provided the U(l) symmetry is 
conserved. Quite often, however, anisotropic interactions 
can change the universality class of the transition and 
open a gap From that point of view, BaCuSi20g 

[9j seems at the moment the most promising candidate 
for the observation of a true BEC quantum critical point 
(QCP) In addition, this system exhibits an un- 

usual dimensionality reduction at the QCP, which was 
attributed to frustration between adjacent planes in the 
nominally body-centered tetragonal structure [ll|]. The 
material also exhibits a weak orthorhombic distortion at 
~90 K which is accompanied by an in-plane IC lattice 
modulation [12| . This structural phase transition affects 
the triplon dispersion, and the possibility of a modula- 
tion of the amplitude of the BEC along the c-axis has 
been speculated based on low held inelastic neutron data 

In order to get a microscopic insight of this system, we 
performed 29 Si and 63 - 65 Cu NMR in BaCuSi 2 6 single 
crystals. Our data in the gapped phase reveal that the 



structural phase transition which occurs around 90 K not 
only introduces an IC distortion within the planes, but 
also leads to the existence of two types of planes alter- 
nating along the c-axis. From one plane to the other, the 
intra-dimer exchange coupling and the energy gap for the 
triplet states differs by 16 %. Exploring the vicinity of 
the QCP in the temperature (T) range 50-720 mK, we 
confirm the linear dependence of Tbec with H — H c \ as 
expected for a 2D BEC. Our main finding is that the av- 
erage boson density n in the BEC is strongly modulated 
along the c-axis in a ratio of the order of 1:5 for every sec- 
ond plane, whereas its local value n(R) is IC modulated 
within each plane. 

NMR measurements have been obtained on ~10 mg 
single crystals of BaCuSi2 06 using a home-made spec- 
trometer and applying an external magnetic field H along 
the c axis. The gapped phase was studied using a su- 
perconducting magnet in the field range 13-15 T and 
the temperature range 3-100 K. The investigation of the 
FIMO phase was conducted in a 20 MW resistive magnet 
at the GHMFL in the field range 22-25 T and the tem- 
perature range 50-720 mK. Except for a few field sweeps 
in the gapped phase, the spectra were obtained at fixed 
fields by sweeping the frequency in regular steps and sum- 
ming the Fourier transforms of the recorded echoes. 

Before discussing the microscopic nature of the QCP, 
let us first consider the NMR data in the gapped phase. 
The system consists of S = 1/2 Cu spin dimers parallel 
to the c axis and arranged (at room temperature) on a 
square lattice in the ab plane. Each Cu dimer is sur- 
rounded by four Si atoms, lying approximately in the 
equatorial plane. For Cu nuclei, the interaction with the 
electronic spins is dominated by the on-site hyperfine in- 
teraction. For 29 Si nuclei both the transferred hyperfine 
interaction through oxygen atoms with a single dimer and 
the direct dipolar interaction are important. According 
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FIG. 1: (Color online) Evolution of the normalized 29 Si NMR 
spectra as a function of T in the gapped phase. Below 90 K 
the line splits into two components, each of them correspond- 
ing to an IC pattern. Inset: T dependence of the 1st moment 
(i.e., the average position) for i) the total spectra (squares) 
and ii) the individual components before they overlap (up 
and down triangles). The solid and dashed lines are fits for 
non-interacting dimers. 
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FIG. 2: (Color online) 63 ' 65 Cu NMR spectra of BaCuSi 2 6 
in the gapped phase, well below the critical field. The T de- 
pendence of the high-frequency "satellite" line clearly reveals 
two different copper sites. From their shifts, the two corre- 
sponding gap values have been determined. Inset: field sweep 
spectrum that reveals the IC nature of the line shape for each 
of the two sites. Shading separates the contribution of the 
65 Cu high-frequency satellite from the rest of the spectrum. 
The analysis of the latter part confirms that the observed line 
shape has a pure magnetic origin. 



to the room temperature structure I4i/acd [14j . there 
should be only one single Cu and two nearly equivalent 
Si sites for NMR when H\\c. As far as 29 Si is concerned, 
one actually observes a single line above 90 K, as can be 
seen in Fig. [TJ However, below 90 K, the line splits into 
two components, each of them corresponding to an IC 
pattern, that is an infinite number of inequivalcnt sites. 
This corresponds to the IC structural phase transition 
discovered by X-ray measurements 12 1. At 3 K, when 



T is much smaller than the gap, the spin polarization is 
zero and one observes again a single unshifted line, at the 
frequency v = v$ = 29 jH defined by the Si gyromagnctic 
ratio 29 j. 

On the 63 - 65 Cu NMR spectra recorded at 3 K and 
13.2 T (Fig. [2]), however, one can distinguish two dif- 
ferent Cu sites, denoted A and B. That is, each of the 
6 lines of Cu spectrum (for 2 copper isotopes x 3 tran- 
sitions of a spin 3/2 nucleus) is split into two, which 
is particularly obvious on the lowest frequency "satel- 
lite" 63 Cu line. The whole spectra can be nicely fitted 
with the following parameters: 63 = 14.85 (14.14) 

MHz, n = 0, and K Z !> B) = 1.80 (1.93) %, where v Q is the 
quadrupolar frequency and rj the asymmetry parameter. 
The K zz is the hyperfine shift, expected to be purely 
orbital since the susceptibility has fully vanished. On in- 
creasing T the highest frequency 65 Cu "satellite" lines of 
sites A and B become well separated and both exhibit 
a line shape typical of an IC modulation of the nuclear 
spin-Hamiltonian. Although the apparent intensities of 
lines A and B look different, they correspond to the same 
number of nuclei after corrections due to different spin- 
spin relaxation rate I/T2. Since the satellite NMR lines 



at 3 K (the lowest temperature) are narrow, the modu- 
lation of is negligible, meaning that the IC lincshapes 
visible at higher temperature are purely magnetic. This 
is confirmed by the analysis of the spectrum shown in the 
inset of Fig. [31 which shows that at 8.9 K the broadening 
of the "central" line is the same as that observed on the 
"satellites". Such a broadening results from a distribu- 
tion of local hyperfine fields: 8h z {IL) = A zz (R)to 2 (R) in 
which ^4(R) is the hyperfine coupling tensor and m z (R) 
the longitudinal magnetization at site R. Since j/q(R) is 
not modulated by the distortion, one expects that the 
modulation of A(R) is negligible too, A(R) =A. This 
means that the NMR lineshape directly reflects the IC 
modulation of m z in the plane. 

Keeping constant the vq parameters obtained at 3 K, 
one can analyze the T dependence of the shift K zz (T) of 
each component a — A or B according to the formula 



K°(T) - K«M = A« z m d z (A a ,H,T)/H, 



(1) 



where m z is the magnetization of a non-interacting 
dimer, m z = g c fj,B/(e^ Aa[ ~ 9atIBH ^ kBT + 1) in the given 
T range, g c = 2.3 [15j, and K zz is determined from the 
average line position, i.e., the first moment. The best fit 
was obtained for Aa(b) = 3.7 (4.3) mcV and A^}^ = 
-16.4 T/^b- We assumed that A^ c = A^ c , but the values 
of A depend only weakly with this quantity. The values 
are slightly higher than those determined by neutron in- 
elastic scattering for Q m i n = [tt, tt] (HI, which is normal 
considering our approximate description. However, the 
ratio Ab/Aa = 1-16 is in excellent agreement with the 
neutron result 1.15. Considering the fact that there is 
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FIG. 3: (Color online) Evolution of the normalized 29 Si spec- 
trum as a function of H at fixed T. The colored spectra 
correspond to the BEC. a) T = 50 mK: Instead of a simple 
splitting of the line as expected for a standard BEC, a com- 
plex pattern appears, typical of an IC distribution of the local 
hyperfine field. Inset: H dependence of the 1 st (Mi, squares) 
and 2 nd (M2, circles) moment of the spectra. Mi is propor- 
tional to m z and Mi to the square of the order parameter, b) 
T = 720 mK: The non zero magnetization outside the BEC 
leads to an IC pattern for fields H < H C \(T), where H c i(T) 
is determined from the H dependence of Mi, as shown in the 
inset. Lower inset: T c is linear in H — H c i, as expected for a 
2D BEC QCP. 



no disorder in the system (as Cu lines at low T are nar- 
row), and that X-rays did not detect any commensurate 
peak corresponding to a doubling of the unit cell in the 
ab plane, our NMR data can only be explained if there 
are two types of planes with different gap values. Look- 
ing back at the 29 Si spectra in Fig. [TJ one also observes 
just below 90 K two well separated components, both of 
them exhibiting an IC pattern. They indeed correspond 
to the two types of planes, as the T dependence of their 
positions can be well fit using values close to Aa and 
Ab determined from Cu NMR (inset to Fig. HJ. This 
means that the 90 K structural phase transition not only 
corresponds to the onset of an IC distortion in the ab 
plane, but also leads simultaneously to an alternation of 
different planes along the c-axis [16j |. with intra-dimcr 
exchange in the ratio Jb/Ja — Ab/Aa = 1.16 



Let us now recall what is expected from a microscopic 
point of view in the vicinity of the QCP corresponding to 
the onset of a homogeneous BEC for coupled dimer sys- 
tems. As soon as a finite density of bosons n is present 
(H > H c i = A m i n /5/LtB), a transverse staggered magne- 
tization m± (_L to H) appears. Its amplitude and direc- 
tion correspond respectively to the amplitude and phase 
of the order parameter. At the same time, the longitu- 
dinal magnetization m z is proportional to the number of 
bosons at a given temperature and field, this latter play- 
ing the role of the chemical potential. Due to the appear- 
ance of a static mj_ , the degeneracy between sites which 
were equivalent outside the condensate will be lifted and 
their corresponding NMR lines will be split into two. To 
be more specific, we consider a pair of Si sites situated 
in the ab plane on opposite sides of a Cu dimer. Outside 
the condensate, and in the absence of the IC modula- 
tion, they should give a single line for if || c. Inside 
the condensate the NMR lines of this pair of Si sites will 
split by ± 29 j\A z ±\m± because their A z ± couplings are 
of opposite sign. Obviously, observing a splitting of lines 
requires the existence of off-diagonal terms in the hyper- 
fine tensor. Such terms are always present due to the 
direct dipole interaction between an electronic and the 
nuclear spin, which can be easily calculated. 

Instead of this expected simple line splitting, the spec- 
tra of Fig. [3^ reveal a quite complex modification of the 
line-shape when entering the condensate. The narrow 
single line, observed at 23.41 T at the frequency 1/0, 
which corresponds to a negligible boson density, sud- 
denly changes into a composite line-shape including a 
narrow and a broad component. The spread-out of 
the broad component increases very quickly with the 
field. The width of the narrow component also in- 
creases, but at a much lower rate. Both peculiar broad- 
cnings are related to the IC modulation of the boson 
density n(R) due to the structural modulation. To be 
more precise, a copper dimer at position R has in to- 
tal 4 Si atoms (denoted by k = 0,1,2,3) situated around 
in a nearly symmetrical square coordination. The ab- 
solute values of the corresponding hyperfine couplings 
will thus be nearly identical, and we will also neglect 
their dependence on R. These 4 Si sites will give rise 
to four NMR lines at the frequencies 29 i/fc(R) = vq + 
fi(R) + f2,jfe(R), where z^i(R) = 29 -fA zz gfi B n(R) and 
f2,fc(R) = 29 lA z ±m±(R) cos((/) — kn/2). Note that is 2 ,k 
only exists when the bosons are condensed, that is when 
there is a transverse magnetization m± pointing in the di- 
rection (p. In a uniform condensate m± is proportional to 
y/n near the QCP, since the mean field behavior is valid 
in both, 2D and 3D. We assume that only the amplitude 
of the order parameter is spatially modulated, and that 
mj_ (R) oc y/n(R). The line shape is the histogram of the 
distribution of 29 ffc(R), convoluted by some broadening 
due to nuclei - nuclei interaction. 

Three quantities can be derived from the analysis of 
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NMR lines at fixed T values and variable H: the average 
boson density n(H, T), the field H cl (T) corresponding to 
to the BEC phase boundary, and the field dependence of 
the BEC order parameter (for T close to zero). The aver- 
age number of bosons n per dimer is directly proportional 
to the first moment Mi (i.e., the average position) of the 
line: M x = f^u - vo)f(y)dv = ^^A^g^H.T), 
where the line shape f(v) is supposed to be normalized. 
The second moment (i.e., the square of the width) of the 
line AI2 = [y — i>o — M\) 2 f(y)Av has two origins: 
the broadening due to the IC distribution of (n(R)-n), 
and that due to the onset of m± cx \J n(R) in the con- 
densate. When increasing H at T ~ 0, the condensation 
occurs as soon as bosons populate the dimer plane. This 
is observed in the inset of Fig. [3^ at T = 50 mK. Both 
Mi (n) and M2 (mx) vary linearly with the field and the 
extrapolation of M2 to zero allows the determination of 
H c i at 50 mK. For higher temperatures a thermal pop- 
ulation of bosons n exists and increases with H before 
entering the BEC phase. As a result both Mi and M2 
increase non-linearly with H, as shown in the upper in- 
set of Fig. However, the increase of M2(H) shows 
two clearly separated regimes and allows the determina- 
tion of H c i (T) as the point where the rate of change of 
M2(H) strongly increases due to the appearance of mj_. 
Applying this criterion to all temperatures, we were able 
to determine the field dependence of Tbec (lower inset of 
Fig. St) and define precisely the QCP at H cl = 23.35 T. 
In agreement with the torque measurements (llj , we find 
a linear field dependence. This is the signature of a 2D 
BEC QCP, where T c cx (H - H cl f with cf> = 2/d and 
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This analysis, however, does not take into account the 
specificity of the line shapes, which are related to the ex- 
istence of two types of planes with different energy gaps. 
A careful examination of the spectra clearly reveals that 
they correspond to the superposition of two lines exhibit- 
ing different field dependence at fixed T value. For sake 
of simplicity, we have made a decomposition only for the 
spectra at 50 mK, as shown in the inset to Fig.UJ Clearly, 
one of the components remains relatively narrow with- 
out any splitting, whereas the other immediately heav- 
ily broadens in some sort of triangular line shape. The 
field dependence of Mi of the two components, shown in 
Fig. 21 reveals that they differ by a factor of 5. This is 
attributed to the difference by a factor of 5 in the cor- 
responding average populations of bosons. If there were 
no hopping of bosons between A and B planes, the B 
planes should be empty for the range of field such that 
Aa < 9^bH < Ab. Although the observed density of 
boson is finite in the B planes, it is strongly reduced, 
giving rise to a strong commensurate modulation of n 
along the c-axis. According to [llj, the hopping along 
the c-axis of bosons in the condensate is forbidden by 
the frustration, and can only occur as a correlated jump 
of a pair. However, this argument does not take into 




FIG. 4: (Color online) Using a simple decomposition of the 
spectra into two components as shown in the inset, we deter- 
mined the 1 st moments of the 29 Si lines corresponding to the 
different types of planes A and B. From the slopes of their 
field dependence, the ratio of the average boson density is 
found equal to ua/ub — 5. 



account the IC modulation of the boson density. 

In conclusion, this NMR study of the 2D weakly cou- 
pled dimers BaCuSi2 0g reveals that the microscopic na- 
ture of the BEC in this system is much more complicated 
than first expected. Two types of planes are clearly evi- 
denced, with different intra-dimcr J couplings and a gap 
ratio of 1.16. Close to the QCP we observed that the 
density of bosons, which is IC modulated within each 
plane, is reduced in every second plane along the c-axis 
by a factor of ~ 5. This provides new constraints for 
the understanding of the quasi-2D character of the BEC 
close to the QCP. 
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